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ABSTRACT 

The brightest Ultra-Luminous X-ray source (ULX), ESO 243-49 HLX-1, with 
a 0.2 - 10 keV X-ray luminosity of up to 10^^ erg s~^, provides the strongest evi- 
dence to date for the existence of intermediate mass black holes. Although small 
scale X-ray spectral variability has already been demonstrated, we have initiated 
a monitoring campaign with the X-ray Telescope onboard the Swift satellite to 
search for luminosity-related spectral changes and to compare its behaviour with 
the better studied stellar mass black holes. In this paper, we report a drop in the 
XRT count rate by a factor of ~ 8 which occurred simultaneously with a hard- 
ening of the X-ray spectrum. A second observation found that the source had 
re-brightened by a factor of ~ 21 which occurred simultaneously with a softening 
of the X-ray spectrum. This may be the first evidence for a transition between 
the low/hard and high/soft states. 

Subject headings: X-rays: individual(HLX-l) — X-rays: binaries — accretion, 
accretion disks 
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Introduction 



The most compelling evidence for the existence of intermediate mass black holes 
(IMBHs) comes from the observation of ULXs, extragalactic X-ray sources with bolometric 
luminosities exceeding 10^^ erg s~^ which are located outside the nucleus of the host galaxy. 
These X-ray luminosities - if they are assumed to be isotropically radiated - are up to 
several orders of mag nitude above the Eddington limit for the maximum mass of stellar 



mass black holes (e.g. [Roberts 20071 ). There is still an open debate about whether any ULXs 
host IMBHs. How IMBHs form and evolve is also a subject of intense debate, but they are 
thought to be associated with events such as the implosion of massive stars formed during 
the very first stages of star formation, the collapse of dense star clusters, and the early 



growth of supermassive black holes lying in the center of galaxies (IMiller fc Colbert 20041 ). 



2XMM JO 11028. 1-460421, referred to hereafter as HLX-1, was discovered serendipi- 
tously by XMM-Newton on 23 November 2004 (hereafter XMMl) in the outs kirts of the 
edge-on spiral galaxy ESQ 243-49, at a redshift of 0.0224 ^Afonso et al. 20051 ). Its 0.2 - 10 
keV unabsorbed X-ray luminosity, assuming isotropic emission, was found to be 1.1 x 10^^ 
ergs s~^: one order of magnitude larger than the previously known brightest ULX (Miniutti 
et al. 2006). A second XMM-Newton observation performed 4 years later (on 28 November 
2008 - hereafter XMM2) revealed that the source spectrum and X-ray luminosity had 
changed. From its highest X-ray luminosity and the conservative assumption that the 
observed value exceeded the Eddington limit by a factor of 10, a lower limit of 500 Mq was 



derived for the b 



ack hole (BH) in H LX-1, thus providing the strongest evidence so far that 



IMBHs do exist ( IFarrell et al. 2003 ). 



The evidence for spectral variability in HLX-1 from the two XMM-Newton observations 
was the motivation for monitoring the source with the Swift X-ray Telescope (XRT; 



Burrows et al. 20051 ). The aim of this program is to study its spectral properties as a 
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function of its X-ray luminosity to search for any of the thre e canonical X-ray states identifi ed 



in Galactic black hole binaries (GBHBs) (see for instance 



Remillard fc McClintock 20061 ) 



Observing transitions between these states could thus pr ovide important information on the 



physical nature of ULXs (e.g. iKajava fc Poutanen 20091 ). Despite its lower effective area 



compared to XMM-Newton and Chandra, the XRT can record enough counts in relatively 
short exposure times (few to tens of kiloseconds) to track any significant luminosity and 
spectral changes. This has been recently demonstrated by Kaaret & Feng (2009) who 
reported on a monitoring campaign with the XRT of three prototype ULXs, i.e. Holmberg 
IX X-1, NGC 5408 X-1 and NGC 4395 X-2. No spectral state changes were observed in 
any of these three sources, despite significant changes in their X-ray luminosities (by up to 
a factor of 9). 

In this paper, we present the results of the first two observations of our on-going 
monitoring campaign of HLX-1 with the Swift XRT. Compared with a previous XRT 
pointing performed one year before, this observation revealed not only a dramatic drop 
in the XRT count rate (by a factor of ~ 8) but also a significant hardening of the X-ray 
spectrum as well as a dramatic re-brightening (by a factor ~ 21) occurring simultaneously 
with a clear softening of the X-ray spectrum. 



2. X-ray observations 

HLX-1 was observed with Swift under the Target-of-Opportunity (ToO) program 
on 3 occasions so far: SI = 2008-10-24 (33.5 ks), S2 = 2009-08-05 (19.2 ks) and S3 = 



2009-08-16 (19.1 ks). A 

tool XRTPIPELINE vO.12.3 



1 the Swift-XKT Photon Counting data were processed using the 



0, 



We used the grade 0-12 events, giving slightly higher effective 



^See [htt^777h 

easarc.gsfc.nasa.gov/docs/swift/analysis/ 1 
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area at higher energies than the grade events, and a 20 pixel (47.2 arcseconds) radius 
circle to extract the source and background spectra using xselect v2.4a. The background 
extraction region chosen close to the source extraction region is the same for the three 
epochs SI, S2 and S3. No XMM-Newton sources are present inside the background 
extraction region. The average background count rate is 6.2 ± 2.2 x 10~^ count s^^. The 
ancillary response files were created using xrtmkarf vO.5.6 and exposure maps generated by 
XRTEXPOMAP vO.2.5. We fitted all the spectra within xspec vl2.5.0a using the response file 
SWXPC0TO1 2S6-20070901V011.RMF, which iucludes an im provement of the soft energy response of 



Godet et al. 2009bl ): which is essential for a source as soft as 



the XRT (IGodet et al. 2n09a : 
HLX-1. 

The XRT monitoring revealed that HLX-1 was highly variable over the past 10 
months with a drop in count rate by a factor of ~ 8 from SI to S2 followed by a dramatic 
re-brightening by a factor of ~ 21 from S2 to S3 (see Tabled]). 

Fig. [1] shows the unfolded spectra for SI (black), S2 (green) and S3 (red). The SI and 
S3 spectra were binned to contain a minimum of 20 counts per channel, and were fitted 
using the minimisation technique. For all spectra, we fixed the absorption column at 
4 X 10^° cm"^, the best constraint on Nh derived by Farrell et al. (2009). All the best-fit 
models and spectral parameters are given in Table [H From Fig. [U it is immediately obvious 
that the source has varied significantly in flux. The results from SI when fitted using 



an absorbed power-law appear to be consistent with those of XMMl ( iFarrell et al. 20091 ) 
which showed a steep power-law, but are inconsistent with those of XMM2 (x^/dof=134/18 
using the XMM2 best fit model) even though the two observations were performed only 
a month apart. The addition of a disk black-body component to the power-law model 
does not significantly improve the fit [Ax^/dof = 3.1/2; which corresponds to a F-test 
probability of 30%). The use of a unique absorbed disk black-body does not give a good fit 
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either {x^/dof = 40.7/16). The S3 spectrum is extremely soft. Within the statistics of the 
XRT data, the S3 spectrum is better fitted {x^/dof = 14.4/21) using a unique absorbed 
soft thermal emission than an absorbed power-law {x^/dof = 55.4/21) with a A'^H-value 
fixed at 4 X 10^° cm~^. The sparse counts in S2 (28 counts in the background-subtracted 
spectrum) prevent us from performing any meaningful fitting analysis. Instead, we prefer 
to use spectral hardness-ratios (see below). Nevertheless, folding through the Swift-XKT 
response kernel the SI, XMM2 and S3 best-fit models multiplied by a constant factor gives 
adequate fits to the S2 data with a constant factor of ~ 8.2 x 10"^ {L ~ 8.1 x 10'*° erg s^"*^) 
for SI, ~ 8.3 X 10^2 _ 5 x 10^° erg s"^) for XMM2 and ~ 3.1 x 10"^ _ 3 4 ^ lo^o 
erg s~^) for S3. The numbers in parentheses are the corresponding unabsorbed 0.2-10 keV 
luminosity. In all cases, the S2 luminosity appears to be much lower than those derived for 
XMMl, XMM2, SI and S3 by at least a factor of 10. For comparison, we plot in Fig. [T]the 
S2 spectrum (green) using a power- law with a fixed value of the photon index F = 2. 

The fitting results are suggestive of a series of spectral changes. To investigate this, we 
compared the total number of counts in three energy bands: 0.3-1 keV, 1-3 keV and 3-10 
keV (see Table [2]). There is a significant drop (by a factor of ~ 16) in the soft 0.3-1 keV 
band between SI and S2 when compared to the 1-3 keV band. This suggests a hardening of 
the spectrum between SI and S2. The source was not detected above 3 keV in either the S2 
or S3 data due to an overall low flux level in S2 and a steep spectral slope in S3. From S2 
to S3, there is a significant increase in counts (by a factor of ~ 44) in the 0.3-1 keV band 
when compared to the other bands, indicating that the X-ray spectrum has softened again 
(see also Table [1]). 
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State transitions in a ULX? 



The detection of the same spectral states as those observed in the GBHBs would 
be extremely valuable in constraining the nature of ULXs. Spectra of other ULXs have 
been either interpreted as the source being in a low/hard state o r in a very high state , 



even if the distinction between the two states is problematic (e.g. ISoria fc Kuncic 20081 ). 
It appears that ULXs in the high/soft or thermal state are rare. Thermal components 
with temperatures in the range kT = 1 - 1.5 kev have been observed in the spectra 
of some ULXs for an emitted lur ninosity Lx < 10^^ erg s~^, placing them at the 
extreme end of stellar mass BHs fiRoberts et al. 20021 ). Winter et al. (2006) claimed 
the identification of several ULXs in the thermal state based on the detection of disk 
black-body (DBB) spectral components. However, in all cases, the DBB contribution was 
not significant with respect to the power-law component. Moreover, no state transition 



as seen in GBHBs has been observed in U 
variability (e.g. 



Gladstone fc R.oberts 2009 



Xs even those showin g large luminosity 



Fridriksson et al. 2008i) . Claims for spectral 



state transitions, different from those seen in GBHBs have already been reported (see e.g. 



Liu et al. 2002 



Soria fc Motch 20041 ). In these papers, the authors claimed evidence for 
a high/hard to low/soft transition. Recently, Isobe et al. (2009) claimed evidence for a 
spectral transition in the ULX, NGC 2403 Source 3 from a slim-disk state dominated by a 
~ 1 keV disk black-body to a very high state dominated by a power-law spectrum. 

When compared to other ULXs, HLX-1 is truly remarkable not only because its huge 
l uminosity enables us to claim that it may harbour an IMBH with a > 500 Mq mass 



( iFarrell et al. 20091 ). but also because its luminosity-spectral variability as observed in 
X-rays (see Section 2) shows compelling evidence for spectral variation on short timescales 
(see Fig. [2]) that are con sistent with what has been observed in GBHBs like GRS 1915-1-105 



( iFender fc Belloni 20041 ). Indeed, HLX-1 was likely to be in the very high state in XMMl 
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& SI, while the source was in the high/soft state in XMM2 with a DBB luminosity 
corresponding to 80% of the total X-ray luminosity. In this case, the transition has occurred 
in a 1 month window. In S2, the data suggest that the source may have been in the 
low/hard state. The re-brightening and spectral softening in S3 suggest that HLX-1 has 
returned to a high/soft state brighter than during XMM2, the transition having occurred 
on a relatively short timescale (< 7 days). 

The low-temperature of the DBB component measured in XMM2 and S3 is in t he same 
range as those measured in some ULXs (~ 0.1 — 0.2 keV; see 



Soria fc Kuncic 2008i). This 



thermal component is often interpreted as direct emission from the accretion disk. The 
DBB normalisation i^(XMM2) ~ 29 and i^"(S3) ~ 14.3 (via the inner radius of the accretion 
disk Rin) could then be used to constrain the BH mass. To do so, we assumed that Rm 
corresponds to the radius of the last stable orbit around a non-rotating BH or a rotating 
BH with a maximum angular momentum (i.e. Rin = 0.5 — 3 Rs ~ 1.5 — 9 km with 
Rs = the Schwarzschild radius where M, c and G are the BH mass, the speed of light 
and the gravitational constant, respectively). We found that the derived masses are in the 
range of IMBH masses with M ~ 5.7 x 10^ - 3.4 x 10"^ M© (cos(6'))"3 from XMM2 and 
M ~ 4 X 10^ — 2.4 X 10'^ Mq (cos(6'))~^ from S3, where 9 is the inclination of the disk with 
respect to the line of sight. 

Within the limited statistics of the XRT spectrum, we have attempted to fit the highest 
quality spectrum (S3) with the slim disk model (Kawaguchi 2003). In order to reduce the 
number of free parameters, we freezed the source distance to 95 Mpc and the viscosity 
parameter to three possible values (0.01, 0.1 and 1). The fits suggest a BH mass larger than 
10^ Mq, thus adding another piece of evidence for the presence of an IMBH in HLX-1. 
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4. Conclusion 

We have presented the first evidence for spectral state transitions in HLX-1, similar 
to those seen in GBHBs, strengthening the case for a black hole in the system. Since 
multi-wavelength observations (e.g. IR, radio) of HLX-1 continue to exclude alternative 
explanations for HLX-1, such as a foreground neutron star, or a background narrow line 
Seyfert 1 galaxy (Webb et al., in preparation), HLX-1 still provides the strongest evidence 
for the existence of IMBHs in the Universe. 

O.G. acknowledges funding from the CNRS and ONES. S.A.F. acknowledges STFC 
funding. We also thank Cole Miller for useful discussions. 

Facilities: XMM-Newton, Swift. 
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Fig. 1.— Swift-XKT PC grade 0-12 unfolded spectra of HLX-1: SI (black) and S3 (red). 
The solid lines correspond to the best-fit models (see Tabled]). For comparison, we plot the 
S2 spectrum (green) using a power-law with a fixed value of the photon index F = 2. We 
used the xspec command setplot rebin 7 7 to improve the visual aspect of the S2 spectrum. 



-13- 



— 1 1 1 r 



i^MMl (VHS) 




T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 



\ S2 (LH) 



_L 



.0 0.5 1.0 1.5 2.0 2.5 

(1-10 keV) / (0.3-1 keV) 



3.0 



Fig. 2. — Hardness intensity diagram using the Swift and XMM-Newton data. The XMMl 
and XMM2 points were computed by convolving with the XRT spectral response the best-fit 
models in XMMl and XMM2, respectively. All the errors are 1 a errors. The S2 0.2-10 keV 
unabsorbed luminosity was derived by folding through the Swift-XBT response kernel the 
SI, XMM2 and S3 best-fit models multiphed by a constant factor (see Section 2) and taking 
the lowest (L ~ 3.4 x 10^° erg s~^) and highest (L ~ 8.1 x 10^*^ erg s~^) values. 
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Table 1: Summary of the X-ray spectral parameters for the best model fits. 



Observation 


Observation" 


Model 


Spectral 






Number 


starting date 




parameters 


(0.2-10 keV) 
(xlO^^ erg s~^) 




SI 


2008-10-24 


ABS PL^ 


r = 3.4 ±0.2 


9.9 ±0.5 


20/16 


S3 


2009-08-16 


ABS DBB'^ 


kT = 0.26 ± 0.02 keV 


11.0 ± 1.0 


14.4/21 



"The SI observation was performed between 2 XMM-Newton observations: XMMl = 2004-11-23 and XMM2 
= 2008-11-28. 

''The unabsorbed 0.2-10 keV luminosity was computed assuming a source distance of 95.5 Mpc and using the 
WMAP cosmology. 

■^ABS PL = absorbed power-law; ABS DBB = absorbed disk black-body. 



djYXMM2 ^ 4 X iq20 ^^-2 ^j^g |-,gg|- constraint on Nh we derived from the XMM2 data (|Farrell et al. 20091 ). 



It is the sum of the Ga lactic column absorption in the direction of HLX-1 N§"-^ = 2 x 10^" cm ^ 
(jPickev fc Lockman 1990( ) and the intrinsic column absorption along the line of sight. 
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Table 2: Comparison of the number of counts in 3 energy bands: 0.3-1 keV, 1-3 keV and 3-10 
keV. The last column gives the source count rate in the 0.3-10 keV band. All the numbers 



derived using 


a 20 pixel radius circle for the source 


and the background. 


Observation 


Exposure 


Counts 


Counts" 


Counts'* 


Count rate 


Number 


time 

(ks) 


(0.3-1 keV) 


(1-3 keV) 


(3-10 keV) 


(0.3-10 keV) 
(cts s"^) 


SI 


33.45 


274 ± 17 


97± 10 


18 ±5 


1.2 X 10^2 




(19.17) 


(157 ± 13)^ 


(56 ±7)^ 


(10 ±4)'' 




S2 


19.17 


10 ±4 


18 ±5 


Q + 2 


1.5 X 10-3 


S3 


19.09 
(19.17) 


437 ±21 
(438 ±21)^ 


185 ± 14 
(186 ± 14)'' 


Q + 2 


3.3 X 10-2 



"The errors quoted above are 1 a error. When the number of counts in a given energy band is less than 
20, the 1-g errors were computed using the following formula: cr = 1 + \/N + 0.75 instead of cr = 
(jsee Gehrels et al. 19861 ). 

''The number of counts are computed using an exposure time of 19.17 ks to facilitate the comparison between 
all the Swift observations. 



